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Abstract 

We derive the relation between the scattering phase shift and the two-particle energy 
in the finite box, which is relevant for extracting the strong phase shifts in lattice QCD. 
We consider elastic scattering of two particles with different mass and with non-zero total 
momentum in the lattice frame. This is a generalization of the Liischer formula, which 
considers zero total momentum, and the generalization of Rummukainen-Gottlieb's formula, 
which considers degenerate particles with non-zero total momentum. We focus on the most 
relevant total momenta in practice, i.e. P = (27r/L) and P = {2ti / L) {e^ + Cy) including 
their multiples and permutations. We find that the P-wave phase shift can be reliably 
extracted from the two-particle energy if the phase shifts for Z > 2 can be neglected, and we 
present the corresponding relations. The reliable extraction of S-wave phase shift is much 
more challenging since Si=o is always accompanied by Si=i in the phase shift relations, and 
we propose strategies for estimating 61=0- We also propose the quark- antiquark and meson- 
meson interpolators that transform according the considered irreducible representations. 



1 Introduction 

The phase shifts for strong elastic scattering of two hadrons encode the basic knowledge on 
the strong interaction between two hadrons, which is non-perturbative in its nature. The 
phase shift 5i is related to the phase between the out-going and the in-going /-wave in the 
region outside the interaction range, and parametrizes our ignorance of the complicated form 
of this interaction. The phase shift indicates whether the interaction is attractive or repulsive, 
what is its strength as well as the range, and it provides the value of the scattering length. 
The knowledge of the phase-shift also serves to determine the masses and the width of the 
resonance that appears in the Z-wave: 61 = 7r/2 at the resonance peak, while the sharpness 



of the rise allows the determination of the resonance width according to Breit-Wigner type 
functional form. In fact, the only feasible method for determining the resonance width on 
the lattice at present goes through the determination of the phase shift. 

Two decades ago, Liischer proposed how do determine phase shifts for elastic scattering 
in a lattice simulation pj. He derived the phase shift relation, that relates the two-particle 
energy ii^ on a lattice of size L and the infinite volume scattering phase shift where 
s = i?^ — and P is the total three-momentum of two particles. If one determines two- 
particle energy E from a lattice simulation with a given momentum P, one can extract the 
phase shift 5{s) at particular s = i?^ — P^ using this relation. Liischer considered only the 
case P = 0. The explicit phase shift relations for higher / at P = ware written down in [2]. 

In order to determine 5{s) at more values of s = E"^ — P^, one better considers also the 
case P 7^ 0. The phase shift relations for this case were derived in [31 El lU |5], but these 
consider only the scattering of two particles with equal mass (mi = m2). 

In this paper we derive the phase shift relations for the general scattering of two particles 
with different mass (mi ^ 1712) and with non-zero total momentum (P 7^ 0). A first step 
in this direction was made by Davoudi and Savage ^ and by Fu [8], where the phase shift 
relation for the irreducible representation Ai was written down: the relation in [7j takes 
into account only the S-wave interaction and neglects all higher partial waves, while the 
relation in [S] takes into account S and P wave. However, the Ai representation is the least 
interesting in practice as it mixes 5* and P-wave phase shifts in one relation if P 7^ and 
mi 7^ m2, making it difficult to reliably extract any of the two. We derive the phase shift 
relations also for the other irreducible representations entering in S or P-wave scattering 
with total momentum P = (27r/L)e2 and P = (2'K/L){ex + e^): these representations do 
not mix S and P-wave phase shifts. We also propose the form of lattice interpolators that 
transform according to these irreducible representations. 

The analogous case of a moving bound state, which is composed of two particles with 
different mass, has been recently explored in [71 [9] . The corresponding finite volume correc- 
tions for the S-wave interaction of two particles has been derived in non-relativistic quantum 
mechanics [H] and in quantum field theory [7|. 

There have been a number of lattice simulations that extracted the phase shift using 
P = 0, or considering P 7^ but mi = m2. The vrvr scattering with 1 = 2 has been studied 
most frequently; see for example [IHl E] for ^-wave at P 7^ and \12\ for D-wave at P = 0. 
The p resonance in vrvr scattering with / = 1 is the only resonance that has been clearly 
observed in the lattice studies of the phase shifts, which allowed the determination of its 
mass and width [131 HH HSl HSl [IZl [18]. The prehminary results for the challenging rni 
scattering with 1 = were presented in [191 120] • The Ktt phase shift was extracted from the 
ground state with P = [21], and the results on / = 3/2 are more reliable than the 7 = 1/2 
ones. The analytic studies of the phase shift relations with P = that may reveal the nature 
of the scalars mesons in these channels were presented in [22]. The pvr scattering and the 
related ai resonance ware simulated in [23] , while the corresponding phase shift-relations for 
the scattering of unstable particles was analytically studied in [21]. The preliminary results 
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for channels including charmed and charmonium states were presented in [23], while D*Di 
scattering relevant for Z^(4430) was simulated in [26]. Recent review of the applications, 
including also baryons, multi-particle interactions and bound states was presented in |27j . 

There are however many interesting channels, where two scattering particles have dif- 
ferent mass, and the simulations at non-zero total momentum would provide the valuable 
information on the corresponding phase shifts. To our knowledge, the phase shifts have not 
been extracted from lattice in such a case, and we provide analytical tools that would enable 
that in the near future. 

In Section 2 we first consider two non-interacting particles in the finite volume, then 
we consider the interacting particles and write down the general phase shift relation. In 
Section 3 we simplify the general phase shift relations by considering the discrete symmetries. 
First we focus on the case of total momentum P = {2Tx/L){ex + e^), write down the phase 
shift relations for three irreducible representations that appear in S or P-wave scattering, 
discuss the strategies for extracting the phase shifts (5i=o,i and provide the quark-antiquark 
and meson-meson interpolators that transform according to these representations. Then we 
repeat the same steps for the total momentum P = {27r /L)ez- We end with conclusions. 
The Appendix provides the derivation of the expression for the generalized zeta function 
Z^{1] q^) for nil 7^ "^2, that is appropriate for numerical evaluation. 



2 Two particles in a finite volume 

We consider a square lattice box of volume with periodic boundary conditions in all three 
spatial directions, while time extent is infinite. We assume continuous space-time and we 
do not consider discretization errors due to the finite lattice spacing a in actual simulations 
with a given action. There are two particles with total three-momentum P in such a box, 
and the total momentum has to satisfy the periodic boundary condition 

27r 

p = Pi + P2 = — d, dez'. (1) 

The main task is to derive the total energy E of the these two particles, where E refers to 
the energy measured by the observer that is at rest with respect to the lattice frame (LF), 
i.e. lattice square box. First we consider the non-interacting case, which is trivial. Then we 
turn to the interacting case, where the energy E depends on the scattering phase shifts 5i 
in the /-th partial wave. This relation will ultimately allow for the determination of 5/ from 
the energies determined by lattice simulations in a finite box. 

The scattering in partial wave I refers to the center-of-momentum frame (CMF), which 
moves with the velocity 

- = l • -' = 71^ 

with respect to the lattice frame. Therefore we need to consider the physical system in CMF, 
where the quantities will be denoted by *. The Lorentz transformation between two systems 
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is performed by 7, which acts on a general vector u as 



7U = 7U|| + U_L , 7 U = 7 U|| + U_L , U|| = PI V , U_L = U - U|| (3) 

so it preserves the component perpendicular to v and modifies the component parallel to v. 
The lattice square box is deformed to some general parallelepiped, and its shape depends 
on the direction of P. The two-particle wave functions in CMF will "see'' the lattice box 
in the shape of this parallelepiped and the technical difficulty is that the periodic boundary 
condition on the CMF wave functions has to be enforced with respect to this parallelepiped. 



2.1 Non- interacting case 

In the non-interacting case there are two major simplifications: the momenta of the individual 
particles also satisfy the periodic boundary condition and the energy is the sum of individual 
energies 

27r 27r 



E^sJ^l + ml + -^pi + mi , pi = — n , p2 = — n' , n, n' G Z'' . (4) 

This already provides the two-particle discrete energy spectrum in absence of interactions. 
We will also need the values of the allowed momenta p* in the CMF frame 

P* = pt = -p; P* = |P*| (5) 
to study the scattering, and we extract it from pi = using 
Pi = 7(P* + ^El) , p2 = 7(-P* + ^E;) so 

p* = 7-^pi - wE{ = 7-^[pi - 7vi?i*] = 7"' - ^ + E*^ ] 

p* = (pi - M p) , (6) 

where in the second step Ei — (£'*/2)[l -|- {m\ — m\)/E*'^] is expressed in terms of energy in 
CMF 

and we have defined coefficient A 

2 2 

mf — 777.9 / X 

^^1 + ^^ (8) 

which is different from 1 only when 7771 7^ r?72. We express the values of p* in terms of the 
dimensionless CMF momentum q 

P* ^ (9) 
and the allowed values of q in the non-interacting case are 

q = r r e Pd (for non — interacting case) (10) 
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where the r G is set of vectors given by the mesh obtained combining ( 6|[9p!0 ) and 



■,-1/ 



1 



^[n - ^ A d) } , ne Z'-^ 

which agrees with [71 |H|^ The symmetries under which this set of points is invariant will 
play a major role later on. The equality q = r (10) will be modified by the two-particle 
interactions in the finite volume. 



'11^ 



2.2 Interacting case 

Now we consider the elastic scattering of two interacting particles with spin in a finite 
box using relativistic quantum mechanics, along the lines of Rummukainen-Gottlieb that 
consider mi = m2 [3], and Fu that presented the analogous derivation for m2 [H]- 

For the case mi = m2, the quantum mechanics result of |3] was subsequently reproduced 
using the Bethe-Salpeter equation [6] and using the quantum field theory [1]. The phase 
shift relations from all approaches agree when one neglects the terms, that are exponentially 
suppressed with the box size L in the quantum field theory. The phase shift relations derived 
here can therefore be applied if L in the simulation is large enough that the terms of the 
order of e~™""^ can be neglected. 

We need to find the two-particle energies E in the finite box in the presence of the 
potential l^(x*), which depends on their relative distance x* = x^ — Xg. The strong potential 
between two hadrons in not known ab-initio in QCD, so one can not analytically calculate 
the eigen-energies E which satisfy Hip{xi,X2) = Eip{xi,X2)-, but rather determines eigen- 
energies E in lattice QCD, which incorporates fundamental QCD interactions. 

However, one can analytically consider the two-particle wave functions in the exterior 
region, where the potential drops to zero 

\/(x*) = for |x*|=x*>i? (12) 

and we assume the interaction is of finite rang^i? < L/2. In the exterior region, the two- 
particle wave function will satisfy Hfree4'free{xi,X2) = Eipfree{xi,X2) with the same eigen- 
energy E as in the case of Hamiltonian H with interactions. The relation Hfree'^free{xi, X2) = 
Eipfree{xi,X2) in the exterior region of CMF has a form of the well-known Helmholtz equation 

{V'+p*')<Pcm{^*) = x*>R (13) 



= + mf + + m| = 7"^ E (14) 

and the total energy E* is just a sum of both individual energies in this region. 

The only effect of the interior region x* < R on the free solutions in the exterior region is 
that 0cA/(x*) will depend on the phase shifts Si{p*), which are related to the phase between 

^We have a different sign than [5] in Pd, but both signs lead to the same (infinite) mesh of points. 
^Presence of the exterior region is not necessary in the quantum field derivation [3]. 
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the out-going and in-going /-wave and parametrize our ignorance of the exact form of the 
potential V{x*). A free solution in the exterior region with momentum p*, which is shifted by 
phase shift Si{p*), will satisfy the periodic boundary condition only for some specific values 
oip*, which fulfill certain relation between p*, 6i{p*) and L. This relation is the analog of the 
Liischer formula we are looking for: it will provide Si{p*) if one determines the momentum 
p* in the exterior region on a lattice of size L. This exterior momentum p* is extracted via 



relation (14) from the energy E of two strongly interacting particles, where E is directly 



measured with a lattice QCD simulation in a box of size L. 

Before imposing the boundary conditions in the finite volume, let us review the familiar 
solutions of the Helmholtz equation for given p* in the infinite volume 



')CM 



y) = Yl yU^, V) W)3i{v*x*) + hi{p*)ni{p*x*)] , x*>R, (15) 



which apply in the exterior region. The phase shift Si{p*) in the continuum is commonly 
defined through the ratio of the out-go in^ /-wave ji — irii and the in-going wave ji + irii with 
momentum p* 

^2is,(p*) ^ ai{p*)+^bi{p*) ^ 
ai{p*) - ihi{p*) 

We will use the same definition of the phase shift in the finite volume, but there the wave 
function will not be so simply expressed in terms of ji and ni due to the boundary conditions. 

Now we turn to the solutions 0c'Af(x*) in the exterior region, that satisfy the Helmholtz 
equation and also the boundary condition at finite L. We consider the case of the periodic 
boundary condition in the lattice frame 

?/'(xi, X2) = i){yii + niL , X2 + n2L) (17) 

which is most commonly used in the actual simulations. These boundary conditions impose 
that 0c'A/(x*) need to satisfy the so-called d-periodic boundary condition, which was derived 
by Fu [8J for P 7^ and mi 7^ 1712 using the Lorentz transformation between the two frame^ 

0cm(x*) = (-1)^""^ 0cm(x* + 7nL) n e . (18) 

A simple example, that satisfies the Helmholtz equation and the d-periodic boundary 
condition, is the Green function 

gik-x* 



G^(x*,p-)=7-^L-3 Y: (19) 



T 27r „ 
fc=-^r, rePd 



•^We apply definition of ni{x) cos{x — Itt/2)/x which agrees with f28^ and T], but differs in sign with 
commonly used definitions. 

^Fu [8] has a different sign here, but this represents exactly the same boundary condition since (—1)"^" '* = 

2 _ 2 

exp(±i7rn • d) exp(zfcz7r '"^,™^ n • d) where each of two signs can be + or — . We choose different sign than 



Fu as it is more in line with our definition of (11) 
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where Pd is a mesh of points (11). Other solutions that satisfy the Helmholtz equation 
and the boundary conditions ( 18 ) are [U El E] 



id/,,* _*2 



(20) 



The solutions Gf^ form a complete basis, as shown by Liischer for d = [Ij, and the general 
solution 0c'm(x*) can be expanded in terms of them 



5cm(X 



X* > R . 



(21) 



Lm 



The technical difficulty arises from the fact that the solutions (21) satisfy boundary 
conditions (18) related to the parallelepiped in CMF, while the phase shifts 6i are related to 
the coefficients in front of the out-going ji — irii and in-going ji + irii spherical waves as in 
the infinite volume (16). So one needs to express the d-periodic solutions (20) in terms of 
the spherical Bessel functions ji and ni. That can be fortunately done in analogous way as 
performed for d = by Liischer [1] and we omit the derivation here 



An 



oo V 



l'=0 m'=-l 



(22) 

where ■MfmVm'i'f') calculable matrices for given /, m, /', m', d, q = Lp*/{2Tc) and A (8) 
and the explicit expression will be given in the next section. The 6 and ip are polar angles 
of X*. 

The general solution in the exterior region (pcAii'^*) (21) is obtained by inserting (22), 
and one needs to relate this to the form ( 15 ) in order to extract the phase shifts defined by 



diel 



?cm(x 



y^vim- "7^- ni 

Z — / 47]- 

Im '- I' ,m' 

X^c/m Yi^{e,^) [ai{p*)ji{p*x*) + bi{p*)ni{p*x*)] , 



X* > R . 



(23) 



By equating the terms in from of YimUi and Yimji we get two relations 



Vim = Clm Ol[p ] 

4:71 



(— i)''(p*)''+^ 



I'm' 



47r 



and Vim can be expressed from the first relation and inserted into the second 

Cl'm' [bl'ip*) Mvm'Mi'f) - 0,l'{P*)^ll'^mm'] = . (25) 

l',m' 

^For our purpose 3^;m(V) was most conveniently applied if both G"^(x*,p*^) and yim{-^, -g^, gpr) are 
expressed in terms of the cartesian coordinates. 
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This hnear system has nontrivial solution for q/^' only if 

det{BM - A) = Aim,Vm' = ai{P*)^ll'^mm' , Bim,l'm' = bi{p*)5w5mm' (26) 

where A and B are defined as diagonal matrices related to coefficients ai and hi |r] and they 



finally provide the information on the phase 5i defined by (16) 



2if A + iB 



By dividing (26) by det{A — iB), which is non-zero [T], one obtains the final relation between 
the diagonal matrix e^**^ and (in general) non-diagonal matrix M 



det[e''\M - z) - (M + z)] = , Ij < Imax (28) 

Mim,l'm' =Mf^i,^,{q'^) , [e^^^]lm,l'm' = 6^*'^'^^ ^5u>5mm' ■ 

This condition is the heart of the phase shift relation and relates the energy E (or q) measured 
on the lattice to the unknown phases 5i{p*) via the calculable matrix elements M.f^i,^i{q^), 
that will be given in the next subsection. The energy level E will provide the information 



on the phase shift 6i{p*) at CMF momentum p* , that is related to E via (14). 



If (5; = for / > Iraaxi the relation (28) needs to be satisfied for the truncated square 
matrices with 1,1' < Imax, as shown in [1]. 

The determinant of the block-diagonal matrix is a product of determinants for separate 
blocks. So the determinant condition will get simplified when M will be written in such basis 
that leads to a block-diagonal form of M and therefore block-diagonal form of e^*^(M— i) — 
(M+i). 

2.3 Definitions of Mf^i,^, and Zf^ for mi m2 

Finally we write down the explicit expression for ■Mfmi'm'il'^)' that were introduced while 



expanding Gf^ in terms of ji and rii ( 22 ) [U, El [8] 

■^'"^/".'(^') - S EiT^'(l'^')^'-WW, (29) 

- (-i)'"'.'-+V(2/ + i)(2j + 1)(2/' + 1) (^ll J^,) (^^^ll 

where Cim,js,i'm' is expressed in terms of the 3j-Wigner symbols. The modified zeta function 
is defined as in [TJ |8] 

Z'™(-i«')-E7ni. (30) 
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where Pd is the mesh of points defined in (11), and is defined in (20). In the special case 
d = 0, the definition of the zeta function ( |30[ ) agrees with the original definition by Liischer 
|lj. The zeta function depends on I, m, = (Lp*/27r)^, d and A ([s]). The Zqq is finite only 
for s > 3/2, but the divergence is not physical as it cancels in the difference between the 
finite and infinite volume result, as explained in the Appendix A, where Zqq is obtained by 
analytical continuation from s > 3/2 to s = 1. The Z^{1; q^) is finite for / ^ 0, but the sum 
(30) converges to slowly for practical evaluation. In Appendix A we derive an expression 



that is suitable for the numerical evaluation and reproduces the known result in the special 
case mi = m2 [31 El [15] . 



2.4 General form of M for /, 



1 



In this paper we are especially interested in extracting the 5* and P-wave phase shifts 
5;=o,i(p*) for the scattering of two particles with different mass. This problem is signifi- 
cantly simplified if the scattering phases for the partial waves with / > Imax = 1 are small 
and can be neglected, i.e. we will assume that 5i>i = 0. This is generally true for small p*, 
where higher partial waves are generally suppressed, and is often true also for moderate p* 
if there is no D-wave resonance in the vicinity. In this case the matrix M is 4 x 4 matrix in 
the basis Im = 00, 10, 11, 1 — 1 and the expression Mfmi'm' (29) leads to the following form 
for general d 



M 



00 
10 

11 

1 - 1 



00 



10 



-iy/3wiQ Woo + 2^20 
\ iy/3wn -a/3u72i 



11 

i\/3wii 

Woo - W20 
— \f^W22 



1-1 

Woo - W2Q I 



(31) 



where we defined wi^n to simplify the notation 

1 

Win 



-'Im 



(32) 



3 Consequences of the discrete symmetries of mesh 



d 



3l| are zero for particular choices of d, which becomes 



Some of the matrix elements M.frn,i>m' 
apparent after explicit numerical evaluation of the wim oc Z^{l;q'^) given in Appendix A. 



This is a consequence of the discrete symmetries of the mesh (11 ) in CMF. It is helpful to 



first study these symmetries and determine the texture of the matrix (31) for particular 



d before inserting M to the master determinant condition (28). 



So we explore in this section the group G of the symmetry elements R that leaves the 

lese symmetries for two 
or permutations. There 



mesh Pd (11) invariant. Later on we study the consequences of these symmetries for two 

and d oc 



most useful types of non-zero momentum d oc and d oc e^: + e^,, 
are several reasons why symmetry consideration will be helpful: 
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It will indicate which are zero, purely real or purely imaginary, as indicated above. 



The Helmholtz solution ( |l9[ ) is invariant under the transformations R E G due to 
the sum over v E Pd- The other solutions Gf^ are generated by applying 3^im(V), so 
they transform like Yim [H E] 



I 

GUR^\p*^) = J2 ^w(^) Gt'{^\f') , (33) 

m'=-l 

where D^^^i{R) is defined as a representation in the bases Yim 

YUR^) = E (34) 

m'=—l 

and is in general reducible. So the solutions Gf^ form a representation of group G, 



which is in general reducible. The energy eigenstates in (21) are certain linear combi- 
nations of Gf^ that transform according to the irreducible representation of the G; the 
representation Dmm'i-R) transformation ^ (33) has the irreducible block-diagonal 
form in this basis. 

The same linear combinations of Yim will also lead to the block-diagonal form of M, as 
shown by Liischer (see section n5.3 of [Ij). We will therefore write down M and search 
for the basis that leads to the block-diagonal form|^ It will turn out that the resulting 
basis indeed corresponds to the irreducible representations of D ( 33p4 ). 



The determinant condition ([28j) is greatly simplified in the bases where M is block- 
diagonal since the determinant of the block-diagonal matrix is a product of deter- 
minants for separate block^ In this case the determinant condition is simplified to 
analogous conditions for separate blocks (i.e. irreducible representations). 

The lattice interpolators, that are written down in the lattice frame, have to transform 
according to the irreducible representation of the group G after transformed to the 
CMF frame. We will provide useful examples of quark-antiquark and meson-meson 
interpolators, that satisfy this property and be used in the actual simulations to extract 
the phase shifts. 



All these reasons prompt us to consider the symmetries of the mesh Pd (11) for separate 
cases of d. 



^For higher Imax it is probably easier to first determine the basis that makes representation D (34) block 
diagonal. 

^Note that e^*''(M — i) — {M + i) is block diagonal when M is block diagonal, since e^"' is diagonal by 



construction (28) 
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b) 




Figure 1: The mesh (11) for d = + Cy and mi 7^ m2 (A 7^ 1) is plotted in (c), while 



(a) and (b) show the steps how to get it. 



3.1 P = {2t: / L){ex + e^) and consequences of symmetries 

We first explain the case of the momentum d = + ey in detail, since it is general enough 
to illustrate the procedure and since the corresponding group has only few elements. All 
the results can be easily generalized to the case d = N{ex + Cy) or any permutation in the 
direction. 

We first need to determine the symmetry transformations R that leave Pd invariant. The 
mesh Prf (11) can be visualized in Fig. [1J3 and is obtained in two steps: 

1. First the cubic mesh r = u. E m. Fig. [l^ is shifted by — ^v4d oc + Cj, and since 
A ^ 1 for mi 7^ m2 ([s]), the origin is not in the center of the unit cell in xy plane (Fig. 
^). The inversion / with respect to the origin is lost as a symmetry at this stage, so 
the corresponding group G will not contain J; this is a major difference with respect to 
degenerate case mi = m2, when the origin is at the center of the unit cell in xy plane 
and / is the symmetry. We will see that this has important consequences, for example 
that sectors with even and odd I will not decouple, which will present challenges in 
certain cases. 



2. In the second step 7 ^ contracts the distances in the direction v oc e^; + and keeps 
the distances perpendicular to that, so the mesh is not modified in z direction. 

The resulting unit cell in Fig. ^ has the form of rhombic prism and the mesh is invariant 
only under four transformations R listed in Table [ij There Id denotes the identity, C„(V) 
denotes a rotation by 2T{/n around V, while o"(V) denotes the reflection with respect to the 
plane perpendicular to V. These four transformations form group C2V1 which has only one- 
dimensional irreducible representations. For the one-dimensional irreducible representation 
"irrep" the transformation i? on a vector u is given by the character x^"^^{R) 



R u 



X 



irrep 



(R) u , x'^^'^iR) = ±1 



for ID irrep . 



(35) 



The characters of the irreducible (^1,25-82,3) and reducible (/m) representations are given in 
Table [T] along with an example of polynomials and vectors u that transform according to 
these representation^ 

®With the change of coordmates — -j^i^x + Sy), e'y ~ -j^iex — &y) and = our notation for 
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respresent. 


dim 


Id 


C2{ex + Cy) 


<y{(ix - ey) 




polynom. 


vector u 


irred. Ai 


1 


1 


1 


1 


1 


I ,x + y 


~t~ 


irred. A2 


1 


1 


1 


-1 


-1 


(/>1) 


(Z>1) 


irred. -B3 


1 


1 


-1 


1 


-1 


z 




irred. B2 


1 


1 


-1 


-1 


1 


x-y 


€x Gy 


Yi=o 


1 


1 


1 


1 


1 


Yoo 






3 


3 


-1 


1 


1 


Yio, Yii, Yi_i 





Table 1: Characters x(-R) = XlflT -^(-^)«« representations D for transformations R G C2i, 
(with principal axis + e^), that leave the mesh Pd in Fig. [T] for d = Cx + ey invariant. 
Representations 2 and S2,3 are irreducible while the representation F'^^ is reducible. 
Example of polynomials and vectors u that transform according to these representations are 
given on the right. 



The functions Yim (34) and also the solutions Gf^ (33) form a representation F^'^ for 
transformations R E C^v, but the 21 + 1 dimensional representation F*^') with m = — /, .., / is in 
general reducible. We will need the number of times (^N^^^^p^ that irreducible representation 
"irrep" enters in F^') [25] 

^ R 

where g is the number of elements R of the group G, while the characters of irreducible 
representations x^^^^^{R) a^id reducible representations x^''KR) are given in Table [l The 
resulting decomposition i^ 

F(°) = A, (37) 

F(i) = © ^3 © ^2 

F(2) = 2Ai © ^2 © ^3 © ^2 

This indicates that the solutions (and also interpolators) that transform according to the 
listed irreps will contain the following partial waves 

53:/ = !, 2,.. 52:/ = l, 2,.. : / = 0, 1, 2,.. A2 : / = 2, .. (38) 

so B3 or B2 in CMF will couple to / = 1, but not to / = 0. These two representations 
therefore provide a rather clean possibility to extract / = 1 phase shift if partial waves with 
/ > 1 can be neglected. The interpolators that transform according to irrep Ai in CMF will 



C2V coincides with the more conventional one, but we stick to our notation as it is more appropriate for 
d = 63, + Cy. Our naming for ^2,3 agrees with [S1[TS] in mi = 7712 hmit, while Bi does not have an analog 
for mi =/= TO2- in [8 is denoted by Bi. 

^The decomposition agrees with but there B3 is denoted by Bi. 
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couple to both / = and / = 1, and there is unfortunately no irreducible representation which 
would couple only to / = 0. This will present a serious challenge for a reliable extraction of 
/ = phase shift in simulations with non-zero total momentum, as will be discussed in more 
detail later on. 



3.1.1 Values of as consequences of the symmetries for d = Cx + & 



y 



he transformations R G 6*2^ leave the mesh P^i invariant, so the sum over r G in Z, 



Im 



(30) can be replaced by the sum over r' = Rx 



Zira{q ) = 2^ _ g2). = 1^ (r'^ - q2)s =2^ (r^ - q2)s smce |r I = |i?r| 

^ r'=Rr ^ 



(39) 



which will have important consequences for some R. Now we can list the properties of Z^ 
for d = + Cy-. 

• Zf _^ = (— 1)™(Z^'^)* which is the consequence of the analogous relation for Yim- 

. Zi^ = t-^{ZfJ* since R = a(e. - e,) G C^.: 

The reflection R = a^e^ — Cy) with respect to the plane perpendicular to — ty 
transforms x — )■ y ^ x, z ^ z or equivalently 9^9, e^"^ oc x + iy ^ y + ix = 
i{x-iy) oc ie-'^. This transforms 1^^ = f{9) (e^^)™ ^ f{9){ie-'^)"' = /(^)i"^e-*™^ = 



i^'iYim)* and then t9h leads to Zf^ = V^iZf^, 



Rewriting Zf^ = Re{ZfJ + i Im{Zf^) this leads to Im{ZfJ = for m = 0,4,.., 
ReiZfJ = Im{ZfJ for m = 1,5,.., ReiZfJ = for m = 2,6,.. and ReiZfJ = 
—Im{Zf^) for I = 3, 7, .., independent of the value of /. 

This property of Z^ holds for any case where the mesh is symmetric under R = 
o"(ea; — Cy), which is true for all d = 0, Cx+Cy, in case of degenerate or non-degenerate 
masses. 

^hn = for / — m = odd since ^"(e^) G 6*2^: 

The reflection R = a{ez) with respect xy plane transforms 9 ^ n — 9 and y?, 
while YUn -9,if) = {-iy'^^%U9, if) so (g leads to Z^ = {-iT'^Zf^. 



• Note that if inversion / would be in G then Zf^ = {—lyz^ or Z^ = for odd 1. 
This would decouple parts of M for even and odd /. In the present case mi ^ m2 and 
d 7^ 0, so / is not element of G and Zf^ is not zero in general for odd /. 

We verified all the above relations to be true also with the explicit numerical evaluation of 
Zf^ using the expression in Appendix A. For example Zf^{mi ^ rri'i) ^ and has equal real 
and imaginary part, while Z{\(mi = 7712) = as required by the symmetries of D2h for the 
mass-degenerate case [5]. 
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3.1.2 The matrix Mf^i,^, for d = + Cy 

The above relations for or equivalently wim (32) simplify the general matrix (31) to 









00 


10 


11 


1 - 1 


00 




/ 


Woo 





i\/3wii 




10 









Woo + 2^20 








11 






-VSwu 





Woo - W2Q 


\/QW22 


1 - 


1 


V 


i\/3wu 





-\/6W22 


Woo - W20 J 



M = Mt,'r 



since iwi-i = —iwli = —wu due to Re{wii) = Im{wu), and W2-2 
Re(w22) = 0, while woo and W20 are real. 
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(40) 
—W22 due to 



3.1.3 Phase shift relations for d = e^; + 



The phase shift relations are obtained from the determinant condition (28), which gets sim- 



plified when M (40) is written in the basis that renders it block-diagonal. 



Extracting P-wave phase shift from irreps or B2 

The part Yio oc z already presents a separate block, which transforms according to -B3: z 
gets multiplied by —1 for reflection with respect to xy plane and rotation around + Cy and 
stays invariant for other two transformations (see Table [T]). The determinant condition for 
this 1x1 block requires det[e^'^^{M^'^ -i) - {M^^+i)] = or equivalently tan((5i) = l/M^^ 
with = _|_ 2^20, so 



d = e^, + , of C2V ■ tan5i(p*' 



TT 



3/2 



7 Q 



(41) 



This is the flnal relation that allows the determination of P-wave phase shift 5i{p*) from 



the energy E of two particles with total momentum P 
interpolators that transform according to P3. 



27r, 



+ e„) when one uses the 



The eigenvectors of the remaining 3x3 matrix M (40) reveal that one simple eigenvector 
is = -7={—iYii + Yi-i) (X X — y which transforms according to B2 (see Table |l|). The 

corresponding eigenvalue is M^^ = woo — W20 — \/6Im{w22) and the determinant condition 
(|28[ for this 1x1 block gives tan(5i) = , or 

71-3/2 ^ q 

d = Cx+Cy , B2 of C2V '■ tanSiip*) = ^ 



V5 



(42) 

This is another relation that allows determination of Si{p*) when using interpolators in irrep 

B2. 

The phase shift relations for irreps P3 (41) and B2 (42) agre^^with the expressions in 
[5| [15] for the case mi = 7712- The remaining representations, discussed bellow, have different 



10 



Note that in OIH] is complex conjugate of our Zf^^^ (30) 
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role in the mi = m2 and mi 7^ m2 cases. 



Problems and strategies for extracting S-wave phase shift from irrep Ai 



The remaining 2x2 matrix M can not be reduced further and spans the space in the 
basis of vectors uf^ = —iYi-i) cc x + y and u^^ = loo oc 1, which are perpendicular to 

the vectors u-^^ = Yio and v^^ above. The vectors x + y and 1 both remain invariant under 
all four R G and belong to Ai irreducible representation (see Table [l]). The remaining 
2x2 block in the basis u^2 contained in 



B 
ab 



^00 

Y ^ 
■J 10 

-iYii+yi- 
72 



^00 

Woo 

-i\/6w 





11 



i\/6wii 

Woo - W20 + \/QIm{w22] 





^10 





Woo + 21^20 





-^yll+yl-l 







Woo - W20 - VQIm{w22) 
(43) 

We kept the other two 1x1 blocks for completeness, so represents the 



with a, 6 = 0, .., 3 

desired block-diagonal form of M in the basis uf\,u^^ ,u^^ and the superscript "B" refers 



to block-diagonal. The determinant condition (28) is equivalent to determinant conditions 



for three separate blocks and two of those were already written in ( 41p2 ). The determinant 
condition for 2 x 2 block leads to the relation 



[ = e^ + Cy , Ai of C2V 
;2*''°(*'*)(Mo^-z)-(M,^ 



00 



^)][e 



(M- -^)- (M- + ^)] 



\ Al- 



io 



(44) 



Among three phase shift relations ( [4l|42|44[ ), the 5*- wave phase shift 5o enters only in (44), 
so let us discuss the problems and possible strategy for extracting 60 from (44). Suppose we 



determine the energy level E using interpolators that transform according to representation 
Ai. The values of E, A, d provide the values of ^g(g^) in (44) at corresponding 

g = p*L/27r (14), so the relation ( [44) ) presents one equation with two unknowns: So{p*) and 
6i{p*). In order to determine So{p*) from (44) 



+ ey , Ai of C2v 



(45) 



tan5o(p*) 



1 - tan(5i(p*) [woo - W20 + y/^ I'rn{w22)] 



Woo - tan5i(p*) [w^q - W00W20 + VQwoo Im{w22) - 12{Re wn 



one would need to know the value of 6i{p*) at given p* (14). The representations B-^ or B2 
allow determination of 6i{p*), but the problem is that they will in general fix 61 at some 
other value of p*, which is related to the energy E measured for the case of representations 
or B2. Since Si{p*) in (44) is needed at p*, Si{p*) can not be simply used if p* ^ p*. 



This indicates that the reliable extraction of S-wave phase shift will be challenging, since 
there is no irreducible representation that would contain only 60 but not 61. We envisage 
several possible strategies to estimate the value of ^o, which might be used in the pioneering 
simulations along these directions, but it is clear that some of these strategies are not rigorous: 
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1. If there exists a region of p* where Si{p*) is known to be neghgible, the condition 
( 44p5 ) recovers the standard form [71 [8] tan(5o) = 1/woo or 



d = e^ + ey , Aiof : tan5o(p*) = L^/ 2^ ^ '^o(P*) (^6) 

and allows the determination of 5o for that p* [8]. That may be for example possible 
at small p*, where higher partial waves are generally suppressed, but p* has to be away 
from any nearby P-wave resonance where 6i{p*) ^ 7r/2 is not small. 

2. If 5i{j)*) is not negligible, and its p* dependence is expected to be mild, then one can 



estimate 6i{p*) needed in (44) from 6i{p*) using the interpolation between p* and p* . 
In this case 5i{p*) has to be determined for several p* using several representations and 
several total momenta P (for example -03^2 at d = e^; + e^). 

3. In the continuum limit, one expects that the energy levels E determined using different 
irreducible representations will agree for a physical state with a given total momentum 
P and given /. In the past this (near) degeneracy across different representations 
served to determine / (or J) of the resulting states. In our case of interest, we expect 
E^-^ = Ej^-^ = Ej^^ for d = + Cy in the continuum limit, so the resulting p* for 
/ = 1 state will be the same for all three irreps. This allows the determination of 



Si{p*) at the desired p* from -82,3; inserting that to the phase shift relation for Ai (44) 
will finally allow the extraction of Sq{p*). In practice, the equality between E or p* 
from different representations will be slightly spoiled by the discretization errors. This 
procedure would still give a relatively reliable estimate of So{p*) if Si modestly depends 
on p*, i.e. if there is no close by narrow P-wave resonance. 

We expect that more a reliable extraction for S'-wave phase shift of two particles with 
mi 7^ 1712 needs to be performed using a simulation with P = 0. There representation 
Ai mixes 6q only with ^4 and even higher partial waves, which can be safely neglected. 
The drawback of sticking to a single total momentum P = is that one needs to perform 
simulations at several lattice sizes L in order to determine So{p*) as several values of p*. 



3.1.4 Quark-antiquark and meson-meson interpolators for d = e^^^ + Cy 

In order to extract the phase shifts, one needs to simulate the two-particle system on a finite 
lattice and determine its energy E in the lattice frame. To create the two-particle state, one 
may use corresponding two-particle interpolator or the quark-antiquark interpolator, that 
couples well to the two-particle state or the resonance that appears in this channel. The 
interpolators are written down in the lattice frame, but they have to transform according to 
the desired irreducible representation after the transformation to the CMF. In this section 
we write down some simple examples of such interpolators with this property, that may be 
used in the lattice simulations. 

For concreteness, our two particle interpolators refer to two pseudoscalar mesons P1P2 
with masses mi and m2, since pseudoscalar mesons are often stable against the strong decay 
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and therefore their scattering is most interesting phenomenologically. In the continuum, the 
scattering state therefore carries = 0"*" for / = and = 1~ for Z = 1 in our examples. 
We write down examples of interpolators for this case, but this can be generalized to the 
scattering of other types of particles. 

All interpolators will be expressed in terms of the currents {i = x, y, z) 

V,{p) ^ J2 e'^'glx)7.g'(x) P(p) = Yl e^''"9lx)75g'(x) S{p) ^ e^^'^(x)g'(x) , (47) 

XXX 

where one can replace the choices of the 7 matrix in the current with any combination of 
7 matrices and covariant derivatives that gives the same transformation property of the 
current. In this subsection we write all the momenta in units of 2tt/L. 

An example of quark-antiquark interpolator for in the lattice frame (LF) and the 
center-of-momentum frame (CMF) is 

(OS)lf = K(e. + ey) {0%)cMF = K(0) . (48) 

The boost in direction P = Cx + Cy does not modify the polarization in z-direction, while 
the momentum in CMF is obviously zero. Interpolator V^(0) in CMF transforms like Cz, so 
according to -B3 representation in Table [1} Examples of interpolators in B2 and Ai are 

{0%)lf = Vxie, + ey) - Vyie, + e,) (0|)cmf = V;(0) - Vy{0) (49) 
iO%)LF = K(e.. + ey) + Vyiex + ey) {0%)cmf = V;(0) + Vy{Q) 
{0%)lf = Siex + ey) {0%)cMF = 3(0) 

where polarizations along ex — ey and ex + ey do not modify their direction when boosting 
along ex + ey More general quark-antiquark interpolators with non-zero momentum are 
constructed in [301 EI] • 

The two-particle interpolators are linear combinations with momentum choices pij and 

P2i 

{0^r^r%)LF = J] Pl(pi,) P2(p2,) (50) 

j=a,b,.. 

i0^^r%)cMF = ^1(PP ^2(-p*) , P* = Pi, = -P2, = T'IPI " ^P] 

j 

j 

such that they transform according to given irrep in CMF, where momenta p^^- = —P2j = P*j 
are given by ([6]). 

An example of interpolators for i?2,3 with this property is 

(OSJ)lf = Pi(d + u)P2(-u) - Pi(d - u)P2(u) (51) 

A I / for P3 1 

d = ex + ey , u=< r (■ u±a 

y ' \ Cx-ey for P2 j 
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b) z 



c) 



* — * 



• • • 



Figure 2: The mesh (11) for d = and m\ ^ 1712 {A 7^ 1) is plotted in (c), while (a) and 
(b) show the steps how to get it. 

and the momenta of Pi in CMF are ^ 



Pla = d + U 

P16 = d - u 



p: = r\pia - UP) = r'(d + u - Ud) = - u)d + u = cd + u 



= cd — u c = 7 ( 
while the interpolator in CMF is 

(0^^ )CMF = Pl(cd + u)P2(-cd - U) - Pi(cd - u)P2(-cd + u) = 0(d, u) 

Due to the following properties it transforms according to ID irreps ^3,2 
C(d, -u) = -C(d, u) or C(d, Ru) = (P) 0{d, u) 



(52) 



(53) 



(54) 



PO(d,u)p-^ = 0(Pd,Pu) = 0(d, x'''"'^(P)u) 



X 



irrep 



(P) 0{d,u) 



since P G C2V leave d = + Cy unaffected, while Pu = x*''^^p(P)u = ±u for u = (P3) or 
u = Cx — Cy (P2), as listed in Table [ij 
An example of interpolator for Ai is 



(O^f^LF = Pl(d)P2(0) {OZ'^'hMF = Pl(cd)P2(-cd) . 



(55) 



This CMF interpolator is unaffected by all P G since Pd = d, so it obviously satisfies 
Ai transformation properties. 

More general hadron-hadron interpolators are considered in 



3.2 P = {2tt /L)ez and consequences of C^^ symmetries 

Equipped with the knowledge how to handle the momentum d = Cx + Cy, which was more 
general, one can easily consider d = 62, which has more symmetry transformations. 

The mesh Pd ( 11 ) in Fig. [2J3 is obtained from n e by the shift —\Aez and the inversion 
is lost at this stage. Then the lengths in direction are contracted due to_7~^. 

The transformations P that leave this mesh invariant are given in Table [2j and they form 
the group C^v It has four ID-irreps 2, Pi,2 and one 2D- irrep where only Ai and E 
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respresent. 


dim 


Id 


C4(e,) 


C2(e,) 




a{ex + Cy) 


polynom. 


vector u 












a(ej,) 


(y{ex - ey) 






irred. Ai 


1 


1 


1 


1 


1 


1 


1 , z 


0, 


irred. E 


2 


2 





-2 








x,y or Yii,yi_i 






1 


1 


1 


1 


1 


1 


^00 






3 


3 


1 


-1 


1 


1 


^10) Yii, Yi-i 





Table 2: Characters for transformations R G 6*2,; (with principal axis e^), that leave the 
mesh Pd for d = Cz in Fig. |2] invariant. In addition to irreps Ai and i?, 6*4^ has also A2 and 
B1.2 but they do not appear for / = 0, 1 so we omit them. Example of simple objects that 
transform according to these representations are given on the right. 



appear for / = 0, 1 of our interest according to (36) [8J 

r(°) = Ai (56) 

r(^) = Ai®E 

r(^) =Ai®Bi®B2®E. 

So the interpolators that transform according to the listed irreps contain the following partial 
waves: 



E:l = l,2,.. : / = 0, 1, 2, .. 5i : / = 2, .. 52:/ = 2, .. 



(57) 



3.2.1 Values of as consequences of the symmetries for d = 

By comparing to the case of d = e^^ + e^, we find that some relations are still valid, some are 
not valid and there are some additional ones due to the additional elements R G 6*4^: 

. Zf_„^ = i-inzfj* since y^_„ = (-1)™!^:;,. 

• Z^ = i^{Zf^* since a{ex — ey) G The consequences regarding specific values of 
m are listed in the corresponding section for d = Cx + Cy. 

• Z;"^ = if m 7^ 0, 4, 8, ... since C^ie^) G C^v- 

Ci{ez) transforms 6^9, ^ + ^vr, so 1^^ = f\e) e*™^ ^ f\e) e*'"('^+^") = 
/(e)(e-/2)-e^-'^, so Zf^ = {iTZt. 

• Zf^ is not zero for / — m = odd in general since ^(e^) ^ 6*4^ (see Fig. [2|. 

• Zf^ is not zero for / = odd in general, since / ^ C4„: 

Therefore the matrix M will not be decomposed into sectors with even and odd /, so 
(5o and ^1 will again mix in some phase-shift relations [8]. 
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We verified all the above relations also with the explicit numerical evaluation of using 
the expression in Appendix A: in particular Zf^^mi 7^ 7712) 7^ as already found in [8], while 
Zfoirrii = 7712) = as required by the symmetries of D^h for mass-degenerate case. 

3.2.2 The matrix Mf^i,^, for d = e. 



The above relations for or wim (32) simplify the general Ai (31) for d = t 



00 
10 

11 

1-1 



00 10 

/ Woo iV^Wio 

-iVSwio Woo + 2^20 



\ 



11 




Woo - W20 







1-1 

\ 




Woo - W20 j 



where all wio are real as indicated in the previous subsection. 



3.2.3 Phase shift relations for d = e. 



(58) 



The phase shift relation are obtained from the determinant condition (28) using the matrix 



M (58). This matrix is already in the block-diagonal form, that can not be reduced further. 



The determinant condition requires that determinant for each of the block is equal to zero. 



Extracting P-wave phase shift from irrep E 

The determinant condition for each of two 1x1 blocks leads to tan(5i) = with 



Woo - W20 



E of G 



tan5i(]9*) 



3/2 



7 Q 



^00 



V5 



(59) 



Basis vectors Yn and Yi-i form two-dimensional irreducible representation E (Table |2]) and 
interpolators that transform according to one of those two will naturally obey the same 
phase shift relation. Note that the more conventional basis vectors x oc Yn — Yi_i and 



y oc Yii + Yi_i will lead to the same matrix M (58) and therefore the same phase shift 



relation (59) applies for them. Our interpolators in E representation will transform like x 
or y. 



The phase shift relation for irrep (59) agree with the expression in [16] for the case 
mi = 1712. 



Extracting S-wave phase shift from irrep Ai 



The 2x2 block of M (58) spans the basis Iqo 1 and Y^ 



under all R E C2V and therefore belong to irrep Ai (Table 21). The determinant condition 



oc z, which are both invariant 



^^We find an additional factor in front of Wiq with respect to [S]. 
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(28) for this 2x2 block requires 

d = e, , of : (60) 

[e^''^^P*\woo - - {woo + tW^'^'^'^Hwoo + 2w2o - - (^oo + 2w2o + ^)] 



where Wij (32) depend on Zf^{l; q^) and d = e^. 



If we know the energy of two particles in irrep Ai on the lattice, the relation (60) 



presents one relation with two unknowns 6o{p*) and which was already noticed in [8j. 

1 - tan 6i{p*) [woo + 2w2o] 



A reliable extracting of 5o(p*) from (60) 



Cz , Ai of Civ '■ tan5o(p*) 



Woo - tan 5i{p*) [wgg + 2woo^^2o - 3w 



lOJ 



is challenging since one needs the value of at the same p* . We proposed several 

strategies for estimating this in the corresponding section on d = e^, + e^^ and the same 
strategies may be used also for d = e^. The only difference is that the ID irreps i?2,3 are 
replaced by the 2D irrep E. 



3.2.4 Quark-antiquark and meson-meson interpolators for d = Cz 

Examples of quark-antiquark interpolators that transform according to E or A\ in CMF, are 



expressed in terms of currents (47) 



(0f )lf = 14(e.) {O'^)cMF = Vu{0) k = x,y (61) 
{OZU = Vziez) iO%)cMF = VziO) 
{OZ)lf = S{ez) {0%)cMF = S{0) . 

The two-particle interpolators O^^^^, which transform according to E or Ai in CMF, 
are constructed in analogous way as for d = Cx + Cy. The LP interpolator (51) and CMF 



interpolator (53) with d = and u = ek=x,y -L d have the form 

iO^''^')LF = Piicz + ek)P2i-ek) - Piiez - ek)P2{ek) , k = x,y (62) 
{C'e^^)cmf = Piiccz + ek)P2{-cez - e^) - PiicCz - efc)P2(-ce^ + e^) , c = 7"^(1 - ^A) . 

The CMF interpolator transforms like ek=x,y, as can be understood from the discussion for 
d = Cx + Cy case, so {0^^)cmf transforms according to irrep E (see Table |2|. 
An example of interpolator in Ai is 

(0S''^)lf = Pi(e.)P2(0) (OAf')cA^^ = Pi{cez)P2{-cez) . (63) 
The interpolator O'^^ has been applied for the study for p resonance [12] where mi = 



m2 = m.„, wmie U^' 



hile O^^^'^ has been written dowE^^ 



but not employed. 



/ e D^fi in the nii = m2 case, so is anti-symmetrized with respect to both particles in [16j . 
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4 Conclusions 



We derived the relations that allow the lattice QCD extraction of the scattering phase shift 
^i{s) from two-particle energy E in the finite box. We consider the scattering of two particles 
with mi 7^ 7712 and with total momentum P 7^ 0. The simulation of the system with P 7^ 
will be important in practice as it will allow the extraction of Si{s) at several values of 
s = E^- P2. 

We find that the P-wave phase shift can be extracted from the irreducible representation 
E of C^v for P = {27i/L)ez (59), or from the irreducible representations i?2,3 of C2V for 
P = (27r/L)(e,. + Cy) ^j^. To be more specific, these relations allow a reliable extraction 
of Si{s) when s is below inelastic threshold, when Si>2{s) can be neglected and when L is 
large enough that powers of e"'"'"'^ can be neglected. If these conditions are not satisfied, 
one needs to generalize the phase shift relations presented here. 

The reliable extraction of the S'-wave phase shift from a simulation with P = will be 
challenging even if the above three conditions are fulfilled. The reason is that Sq{s) appears 
together with Si{s) in the Ai representation when P 7^ and mi 7^ m2. This mixing happens 
since the inversion is not the symmetry of the two-particle system in CMF. We propose 
several strategies that allow an estimate of 6o{s) at P 7^ in spite of this problem. We 
expect that a more reliable extraction of S'-wave phase shift for two particles with mi 7^ m2 
needs to be performed using a simulation with P = at several values of the lattice size L; 
in this case 60 mixes only with 5j>4 and these can be safely neglected. 

Besides the phase shift relations, we wrote down also the quark-antiquark and meson- 
meson interpolators that transform according to the considered irreducible representations. 
These can be used in actual simulations. 
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Appendix A: Evaluation of Zf^[\\q^) for mi ^ 1712 

In this appendix we derive a form of the generahzed function Z^{1; q^) that is appropriate 
for numerical evaluation. We consider the most general case mi ^ ma, d = ^ and 
general / and m, which has not been considered before. Some parts of our derivation are 
similar to Appendix A of [33], done for / = m = and mi = m2, and to Appendix B of [8], 
done for Z = 1 and m = 0. 

The Zf^{s; q^) for the general case of mi ^ m2 and d = 7^ is defined in (30) 

where the relations for the phase shift depend on Zf^{l]q^) evaluated at s = 1. Here q^ is 
real and can be positive or negative (14). The sum goes over the mesh defined in (11) 
and plotted in Fig. [Tj:; for d = Ca, + Cy and in [2|:; for d = e^. 

The sum is finite at s = 1 for every / and m except for / = m = 0, and we will derive 
the expression that converges faster than (64), and is appropriate for numerical evaluation. 
We will show that sum converges only for s > 3/2 (but not s = 1) in case of / = m = 0. 
The divergence that appears for s = 1 will be exactly equal to the divergence that appears 
in the infinite volume. Since the phase-shift relations depend on the finite volume shift with 
respect to the infinite volume, we will get rid of the divergence by the analytic continuation 
from s > 3/2 to s = 1. 

First we express l/(r^ — g^)* using the definition of the Gamma function Jq°° dt t^~^ e~^"- = 
r(s)/a* and then split the integral to two parts 

1 p1 poo 

= ^ V yiM { / dt f-'e-'^^'-"'^ + / dt t^-^e-^^^^-"^) } . (65) 

r(s) Wo A 

The integral in the second term is finite at s = 1, it is easily evaluated, and renders faster 
convergence than the original sum 



second term = V3^,„(r)— / dtf-^e''^'"-'^"^ ^ y^yin.{r)—, ^. (66) 



The first term (65 ) contains the sum ^rePd -^('^)' which is equivalent to the sum Xlne^a -^('^( 
and we express it using the Poisson summation formula 



/(n)= /rf'^^/We"-"" (67) 
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leading to 



first term = -L [\t f^'e'"" V /n , fn= [ d^x 3^i„(r(x)) e'^l'-WI' + '^^-x ^gg^ 



with r(x) = 7~^(x — ^Ad) (11). We change the integration variable from x to r using 
(fx = det{J)cPr = •yd^r and separate terms that depend only on r using yim{r) = r^Yim{0, (p). 
Applying x = 7r + \Ad (11 ) the term dependent on A factorizes 



/n = 7 e' 



xAnd 



r'^dr e '""^ 







sin 9d9 







2tt 



Yi 



lm\ 



ikr 



(69) 







with k = — 27r7"^n. We insert the well known relation for e" 

oo I' 



jkr 



(70) 



i'=0 m'=-l' 



The integral sin 6 d6 J^^ d(j)YimiO,(j))Yif^,{9,(f)) = Su/dmm' simplifies (69) to 



/„ = 7 47r i-tY (-l)^-d YimiOk, <Pk) / dr e"*''^ r' jz(fcr) . 



,2 „-ir2 „i 



(71) 



The remaining integral can be evaluated with Mathematica 



I I 1 \Aa.A 



U = l{-^y (-1) 



2t 



Y, 



lm\yk-, (Pk) 



3/2 



(72) 



and we apply (fc/2t)' Yim{9k,4>k) = 3^«m(k/2t) = 3^im(-7r7n/t). Inserting this /„ to (|68| we 
get 



first term 



yi 



TT'-fn I 71 



im \ 



t '\t 



3/2 



-(7r7n)2/t 



(73) 



In the case of s = 1, this integral over t is finite for all n except for n = 0. The n = 
divergence occurs only for / = m = since 3^zm(n = 0) oc SioSmo- The term with n = is 
the infinite volume /n=o = / d^xf(x.) analog of X]n/('^) Poisson's formula (67) and 

is finite only for s > 3/2. In order to get rid of the divergence, that cancels in the difference 
between the finite and infinite volume result anyway, we split the n = term in two parts 



1 



T(s] 



dt t'-^/^e"'' 



1 



T(s] 



f dt f'^/^e'"" - 1) + dt 
Jo Jo 



(74) 



s>3/2 1 s^l 



The first integral is finite for s = 1, while the second integral t^ ^^'^dt — g_^/2 
is finite only for s > 3/2, but we analytically continue it to s = 1. 
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Collecting (66) as well as convergent and divergent piece of (68) to get (65), we get finally 



1 /-^\ 3/2 



2 , X / TT 



■\ ' 1 



+ 7 / dt {e'' -j -j=Sio5mo - 7'^Sio5mo 



+ (75) 

re Pa ^ 

which is used for our numerical evaluation and converges rapidly for I, m, d of our interest. 
It is applicable for > and < 0. We verified numerically that this respects all 
the relations listed in the main text, that follow from discrete symmetries at d = + or 
d = e^. 

In the special case mi = m2, our result agrees with the result in [5], which was presented 



for mi = m2 without derivation We also verified that such Z^*^ numerically agrees with 



obtained for mi = m2 via as proposed by 
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